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We perform variational Monte-Carlo calculations to show that bosons in a rotating optical lattice 
will form analogs of fractional quantum Hall states when the tunneling is sufficiently weak compared 
to the interactions and the deviation of density from an integer is commensurate with the effective 
magnetic field. We compare the energies of superfluid and correlated states to one-another and 
to the energies found in full configuration-interaction calculations on small systems. We look at 
overlaps between our variational states and the exact ground-state, characterizing the ways in which 
fractional quantum Hall effect correlations manifest themselves near the Mott insulating state. We 
explore the experimental signatures of these states. 
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I. INTRODUCTION 

We consider the interplay between three paradigmatic 
quantum states of bosons in rotating lattices: Mott in- 
sulators, superfluids, and fractional quantum Hall states. 
The Mott insulator is found when there is an integer 
number of particles per lattice site, and the tunneling 
is sufficiently suppressed relative to the interactions. It 
is an incompressible state, where interactions freeze the 
particles in place. In the standard cartoon, when the 
density of such a system is tuned away from commen- 
surability, the excess particles (or holes) "skate" across 
the frozen Mott sea, forming a superfluid. If the sys- 
tem is rotating, one expects that the collective motion 
of this superfluid will produce a vortex lattice. In 2007, 
Umucahlar and Oktel [l[ argued that when the rotation 
rate is high enough that the number of vortices is com- 
parable to the number of excess particles, then the su- 
perfluid will be unstable to forming a correlated state 
of matter with particles bound to vortices - a situation 
analogous to that found in the fractional quantum Hall 
state. They supported this argument by estimating the 
energy of the superfluid and the correlated state. Here we 
confirm this scenario through more rigorous calculations. 
By using Monte- Carlo techniques we compare the energy 
of variational states describing fractional quantum Hall 
states and superfluid vortex lattices with each-other. We 
also compare these energies with exact results calculated 
for small numbers of particles. We find that there is 
a range of parameters for which the fractional quantum 
Hall states are more favorable than superfluid states. We 
note, however, that the energy differences between these 
states scales as the tunclling energy, and can be quite 
small. 

Most previous studies of analogs of fractional quan- 
tum Hall states in optical lattices have focussed on the 
low density limit, where there are much fewer than one 
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particle per site. In the context of cold atoms, Hafezi, 
S0rensen, Demler and Lukin Q gave an excellent review 
of the basic physics of this limit (including symmetry and 
topology arguments), and argued that one can continu- 
ously deform a Mott insulating state into a fractional 
quantum Hall state by varyingthe strength of an ad- 
ditional superlattice potential Q. They also proposed 
using Bragg spectroscopy to probe these states. Palmer, 
Klein, and Jaksch [|| performed a number of calcula- 
tions focussed on the role of the trap, detection schemes, 
and on inhomogeneities which can spontaneously ap- 
pear in these systems. Bhat, et al. @ carried out full 
configuration-interaction calculations for a small number 
of particles in a rotating lattice with hard-wall boundary 
conditions. Moller and Cooper analyzed the relevance of 
composite fermion wavefunctions to describing these sys- 
tems @ . Nigel Cooper recently produced a review of the 
physics of rotating cold atom clouds including analogs of 
the quantum Hall effect in lattices Q. These, and our 
present study, build on initial works motivated by solid 
state systems ||. 

Translating these arguments to higher densities is not 
completely trivial. The superfluid near the Mott state 
is more complicated than the standard cartoon suggests. 
For example, the mean- field description treats it as a two- 
component plasma of particles and holes, with a small 
imbalance between the density of particles and holes. 
Despite these complications, we find that when the de- 
viation of the particle density from an integer value is 
commensurate with the magnetic flux one can indeed see 
analogs of the fractional quantum Hall effect. 

We start our analysis with the well-known Bose- 
Hubbard Hamiltonian in an effective magnetic field 

H = -t^a\ aj e iA » + - ^fWte - 1), (1) 

where (aj) is the bosonic annihilation (creation) oper- 
ator at site i and hi — a\a,i is the number operator. The 
tunneling is parameterized by t and on-site interactions 
by U. We use the Landau gauge A = (—By,0), so the 
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phases = exp(ie/hc f r * A-dl) acquired when hopping 
in ± x-direction are =F2irai y , where i y is the y coordinate 
scaled by lattice constant a, and in y-direction Aij = 0. 
Here, a = Ba 2 /(hc/e) = p/q is the flux quantum per 
plaquette and we take p and q to be relatively prime in- 
tegers. The single particle spectrum for this problem is 
the famous Hofstadter butterfly [t| . The phase boundary 
between the Mott insulator and superfluid carries signa- 
tures of this single-particle physics [J, [l0l - (l2l | . Away from 
the tips of the Mott lobes, the physics of the superfluid- 
Mott transition of the non-rotating system is in the uni- 
versality class of the dilute Bose gas. Thus we expect 
that phenomena which can be seen in the dilute Bose 
gas will occur there, including the analogs of fractional 
quantum Hall physics which we are exploring here. 

II. VARIATIONAL WAVEFUNCTION 

A. Laughlin State 

We consider the variational ansatz 

l*)= J2 i>(zi,-,z N )al i ...al N \y MI ), (2) 

Zl,...,ZN 

I 



where |*mj) = JJjia])" / Vn^-Wac) is the Mott insula- 
tor state with no par ticles per site and if) is the Laugh- 
lin wavef unction [13| with filling v = 1/m. To de- 
scribe bosons, m must be even. The complex coordinate 
Zi = Xi + iyi specifies the location of the ith particle, with 
i running from 1 to N, where N is the number of excess 
particles. The sum over Zi is a sum over all lattice sites. 
To describe a state with excess holes, we replace a) with 
n. 

To minimize the role of boundaries, the model in ([T]) is 
typically either solved on a sphere or a torus [1, [3] • We 
will work in an L x L torus geometry, corresponding to 
quasiperiodic boundary conditions 



tp(..., z k + L, ...) = tp(...,z k ,...) (3) 
i>(...,z k + iL,...) = e - * 3 "?*"*^...,**,...). 



For these boundary conditions the Laughlin wavefunction 
can explicitly be written as [HJ 



m N 

V^i, zn) = Ne iK * *'e-*. ^ * e -*& E< »? J| ^ [ {Z _ Zp) l] J| [ {Zt _ Zj) l] }™. (4 ) 

13=1 i<j 



Here, AT is the normalization factor, Z = ^ i Zi is N 
times the center-of-mass coordinate, Zp = Xp + iYp are 
the a-priori arbitrary locations of the center-of-mass ze- 
ros. To satisfy the boundary conditions, one requires 
X/3 = ri\L (n\ G Z), K x = 2irn2/L (ri2 € Z), and 
K y = —2,Ti'^2 l pYplL 2 . The quasi-periodic Jacobi theta 
functions are defined by 

oo 

l?l(z,e" rT ) = ^(-l)"- 1 /2 e iTT(n+l/2) 2 e (2n+l)i 2: _ 

— oo 



For our square geometry r = i. This function is odd 
with respect to z and has the following quasi-periodicity 
properties: $i(z + n) = — and "&\{z + ttt) = 
—e~ % ^ T e~ 2%7 -'&\(z). The relation between the flux quan- 
tum per plaquette a — N^/L 2 , filling fraction v — N/N^, 
and excess particle density e = N/L 2 is succinctly given 
by cat — e, where denotes the number of flux quanta 
in the L x L lattice we consider. In what follows, we will 
restrict ourselves to the v = 1/2 Laughlin state (m = 2), 
so that the commensurability requirement between the 
magnetic flux and particle density becomes a = 2s. 



B. Superfluid State 

We will compare the Laughlin state introduced in 
Sec. Ill Al with a Gutzwiller mean field state 

i*^>=n(E ( 5 ) 

where are variational parameters. This wavefunction 
is commonly used to describe the superfluid in the Bose- 
Hubbard model [l5j . It is exact in the non-interacting 
limit and captures the effect of number squeezing. Its 
main deficit is that it does not capture any of the short- 
range correlations in the superfluid. Regardless, the en- 
ergies it produces are good estimates of the superfluid 
energy. In the non-rotating case, the superfluid is trans- 
lationally invariant, and the coefficients P n are indepen- 
dent of i. In our case, where the lattice is rotating, a 
vortex lattice forms, breaking translational invariance. 

Near the Mott lobe, the site occupations are dominated 
by n = tiq and n = no ± 1: that is it is extremely un- 
likely to have more than one extra particle or hole on a 
given site. We therefore truncate our basis to only these 
three values of n. This will also facilitate direct compari- 
son with configuration-interaction calculations using the 
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same truncated basis. We work in an L x L lattice, using 
the boundary conditions which are equivalent to those in 
Eq. ©. 

The numerical techniques for optimizing the /* are well 
documented [Io| , and we will not repeat the detailed dis- 
cussion here. These can be described in terms of a varia- 
tional calculation where one minimizes (\& mfIHoI^ mf) 
with the constraints that the total number of particles M 
and normalization {^mf\^mf) are fixed: this involves 
introducing the chemical potential /j, and a number of 
other Lagrange multipliers. In practice it is more con- 
venient to write H = Hq — (J,M, and follow an iterative 
procedure based upon mean field theory. These two ap- 
proaches are completely equivalent. In comparing ener- 
gies with our other variational state, one must be cau- 
tious and be sure to use (H ) = (H) + fj,M. 



III. EXACT RESULTS ON SMALL SYSTEMS 

A. Approach and Results 

For small systems we can exactly diagonalize the 
Hamiltonian in Eq. (fTJ), taking a configuration- 
interaction approach where we truncate the allowed num- 
ber of particles on a given site to be no , n-o — 1 , or no + 1 . 
For definiteness we take n = 1 : changing this value just 
scales the hopping matrix elements t. For these small 
system sizes we can also directly calculate (^\Ho\^). In 
Sec. IIVI we will discuss larger systems where we need to 
resort to a Monte-Carlo algorithm for calculating this en- 
ergy. 

We consider 12 particles in a 3 x 3 lattice, so that the 
excess particle density is 1/3. We take i/ = 1/2 and 
accordingly the number of quanta of flux per plaquette 
is a = 2/3. Fig. [1] displays the energies (measured in 
units of U) of the first few hundred exact energy eigen- 
states together with the energies of our two variational 
wavefunctions: Eqs. ^ and ([5]). We emphasize that our 
ansatz for the fractional quantum Hall state is not just 
the Laughlin state, where flux is bound to each particle, 
but is rather the coexistence of a Mott state and a Laugh- 
lin state, with flux bound only to the excess particles. 

In Fig. Q] we also show the estimate from Ref. [l|, 
which is supposed to describe the correlated state near 
the Mott insulator, 



AE = Un e - t(n + l)/(a)e, 



(0) 



where the first term represents the on-site interaction of 
excess particles with the Mott insulator and the second 
term is the hopping energy of particles in the Hofstadtcr 
ground state denoted by —tf(a), f(a) > being the di- 
mensionless maximum eigenvalue of the Hofstadter spec- 
trum. Note that t is enhanced by a factor of (no + 1) 
due to the Mott background. No interaction energy is 
included, as it is expected that in this regime the excess 
atoms avoid one- another. It is remarkable how closely 




0.15 



t(U) 



FIG. 1: (Color online) Exact many-body spectrum for 12 
particles in a 3 x 3 lattice with a = 2/3 considering only 0,1, 
and 2 atoms per site (for v = 1/2, excess particle density is 
e — au — 1/3). Also shown by the black solid line is our vari- 
ational estimate of the energy of a fractional quantum Hall 
state of excess particles in the presence of a Mott background. 
Dash-dotted blue line shows the Gutzwiller mean-field super- 
fluid energy for the same density (1+e), corresponding to a 
vortex lattice where the cores are filled with Mott insulator. 
The dashed red line is the estimate of the ground state en- 
ergy from Eq. (BJ, first introduced in p[. For low enough t the 
variational energy of the correlated state of excess particles is 
lower than the superfluid energy. 



this estimate matches the results of the exact diagonal- 
ization for small t. 

For t < 0.13 the energy of our candidate fractional 
quantum Hall state (with optimized Zp) is lower than 
that of the superfluid, while the opposite holds for larger 
t. Our physical picture of this is that as t grows the 
Mott insulator melts, and the density of mobile atoms is 
no longer commensurate with the magnetic field. 

For very small i, the fractional quantum Hall state's 
energy agrees very well with the exact ground state en- 
ergy: this is shown more clearly in Fig.[2jb). In Fig.[2]Ja) 
we show the overlap between our variational state and the 
exact ground state. At low t the overlap is greater than 
95%, but it falls off with increasing t: presumably due 
to the increasing importance of particle-hole excitations. 
The overlap between the ground state and the mean-field 
superfluid (inset of Fig. 0a)) is never large, and their 
energies in Fig. Q] never approach one-another. We be- 
lieve this is due in part to the fact that the mean-field 
state breaks translational invariance, and consequently 



involves a superposition of many eigenstates [16[. A 
quantum superposition of vortex lattices, may in fact be 
a good alternative description of the fractional quantum 
Hall state. 

Given the small difference between the energies of our 
two variational states, one must be somewhat cautious 
about ascribing too much significance to the crossing at 
t ~ 0.13. One also might be concerned that at that value 
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FIG. 2: (Color online) a) Overlap between the u = 1/2 FQH 
+ MI state [|^} from Eq. ([2J] and the exact ground state 
[IGS 1 ), determined from diagonalizing Eq. JTJ) in a truncated 
basis] as a function of tunneling strength t, using the same 
parameters as Fig. [1] Also shown in the inset is the overlap 
between a superfluid vortex lattice and \GS). b) Comparison 
of the variational and exact energies - from Fig. [1] 



of t, both variational states have an energy which is sig- 
nificantly higher than that of the ground state, suggest- 
ing that neither may be particularly good descriptions 
of the true ground state. A third concern is that there 
is no sign of a phase transition in Fig. |2fa): the overlap 
between the fractional quantum Hall state and the exact 
ground state remains above 75% out to t ~ 0.15. Despite 
these caveats, the large overlaps at small t is convincing 
evidence that the ground state at low t is a fractional 
quantum Hall state of excess particles, and it would be 
surprising if the system formed a correlated state at large 
t. 



B. Variational Parameters 



In Fig. [3] we show how the energy of the variational 
state depends on the parameters Zp, which represent 
where "vortices" can be found around which the cen- 
ter of mass flows. The boundary conditions in Eq. ^ 
force the wavefunction to have m = \jv of these zeros 
(in the present case m= 2). In the absence of the lattice, 
the energy is invariant under changing these parameters, 
leading to an m-fold degeneracy of the ground state [13] ■ 
Here this symmetry is absent and the energy depends on 
Zp . Not surprisingly, the overlap between the variational 
state and the exact ground-state is directly correlated 
with the energy. This overlap has a maximum when the 
variational energy has a minimum. 



FIG. 3: (Color online) Variational energy (a) and the overlap 
with the exact ground state (b) as a function of center of mass 
zeros Z\ = Xi + iY\ , Z2 = L — Z\, measured in units of the 
lattice constant. As with Fig.[TJ we consider anixL cell with 
L = 3, flux per plaquette a = 2/3, filling factor v = 1/2 and 
total particle number M = 12. We take K x = 0, K y = 0, and 
t = 0.01(7. The lower the variational energy, the higher the 
overlap. The lowest energy occurs for X\ =Yi = L/2 where 
the overlap is 96.4%. At this point, the variational energy is 
0.318617, which is very close to the exact ground state energy 
of 0.3176(7. 



IV. VARIATIONAL MONTE-CARLO 

Unfortunately the maximum size of the system which 
can be treated by the techniques of Sec. IHIl is quite lim- 
ited. Our preceding results for small system size pre- 
dominantly serve as a guide for physical intuition, and 
cannot quantitatively describe the physics of the infinite 
system. Here we introduce a variational Monte-Carlo 
(VMC) algorithm [l8| in order to calculate the energy 
= (*|.Ho - fJ-M\^>), where M is the total num- 
ber of particles. This will allow us to make a more solid 
comparison of the energies of the superfluid and corre- 
lated states, and draw the phase diagram in Fig. 2) This 
phase diagram illustrates the regions oft — fj, plane where 
either the superfluid or correlated state has a lower en- 
ergy. 

We begin by introducing a basis \R = {zi, • • • , zjv}) 
where the N excess particles are at sites z\ through zn- 
This allows us to write 



RR> 



(R\H\R')(R'\y) 



(7) 



E 



E 



We use a Metropolis algorithm to sample the sum over 
R. Starting from some configuration R we generate 
a new one i?i by attempting to move a single parti- 
cle by one site. We accept the move with probability 
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TABLE I: Results of our variational Monte-Carlo calculation. 
a is the number of flux per plaqette, e is the density of excess 
particles, L is the system size, N is the number of excess par- 
ticles, — (1 + no)tK is the hopping energy per site, where t is 
the hopping matrix element and no is the number of particles 
per site in the underlying Mott state. K is dimensionless. Our 
estimates of the statistical error in K from a binning analysis 
of 80000 samples are given by SK. 



mm{l, Prj/Pro}: we then continue the procedure to 

generate R2, R3, In the resulting Markov chain each 

configuration R will appear with probability Pr. After S 
steps, the energy is then estimated as E$ = Xh=i Eiu/S- 
As is usual, we discard the first few thousand steps so as 
not to bias the sum by our choice of initial configuration. 
We use a binning analysis to estimate the statistical error 
on our sum [19(. 

For each R, we calculate En directly. The Hamiltonian 
only connects a finite number of different configurations 
(those which differ by moving one particle by one site), 
and the sum is straightforward numerically. 

As a further simplification we note that E{fi, t) = 
E (fj,) — (l + no)tK, where E = Un e + U(n Q — l)n /2 - 
fj,(rio + e) is the expectation value of the on-site terms in 
H and — (1 + no)tK is the expectation value of the hop- 
ping energy. K is independent of no, as the only role of 
the Mott background is to provide a Bose-enhancement 
term of (1 + no). By using the Monte-Carlo algorithm 
to calculate K, rather than E, we produce E(p, t) for all 
no, A*, and t at once. 

Table U lists the parameters for which we have per- 
formed VMC calculations. For the smallest system sizes 
(L = 3, 4, 5) we find that the VMC agrees with the di- 
rect calculation of the variational energy. From the chart, 
we conclude that finite size effects are significant in the 
L = 3 cases, but for larger L the differences between the 
energies of the two systems are within a few percent. We 
have not extrapolated to L = 00. 

Fig. H] illustrates our results for a = 1/4. Near the 
constant density line n = 1 + e, there is a region where 
our variational wavefunction has a lower energy than the 
Gutzwiller mean field vortex lattice. This corresponds to 
an incompressible v — 1/2 bosonic Laughlin state above 
the no = 1 Mott insulator. The same argument can be 
advanced for holes by just changing the creation oper- 
ators in Eq. ([2]) with annihilation operators leading to 
a coexisting Mott insulator and FQH state of holes near 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 
t(U) 



FIG. 4: Phase diagram for a = 1/4 and v — 1/2. Bound- 
ary between Mott insulator (MI) and superfluid (SF) states 
is found from a mean-field calculation. Excess particle (or 
hole) density is e — ca> — 0.125. Boundary of the coexistent 
v = 1/2 FQH state of excess particles (holes) and no = 1 
MI state centered around the 1.125 (0.875) constant density 
line is determined from a comparison of VMC and mean-field 
energies. We consider 8 particles in an 8 x 8 lattice in the 
VMC calculation. 

the 1 — £ line, although it is less visible than in the particle 
case. 



V. CREATION AND OBSERVATION 

Several labs currently have the technology to create 
a rotating optical lattice [2(| Hl| , which can be directly 
used to create the system described here. Those exper- 
iments still are far from the Mott regime, but they are 
progressing rapidly. An alternative approach to imple- 
menting Eq. ([TJ is to use a non-rotating lattice, and gen- 
erate the phases on the hopping matrix elements by some 
other means (22j. The most advanced demonstration of 
this technique was from Lin et al. [23j]. 

One of the more promissing schemes for experimen- 
tally observing the incompressible states described here is 
through in-situ imaging of the density profile of a trapped 
gas 0,[2J]. The fractional quantum Hall states should ap- 
pear as extra steps in the density profile near the Mott 
insulator plateaus. Moreover, as the magnitude of the ef- 
fective magnetic field increases, these steps move in pre- 
dictable ways: the density is set by the magnetic flux, 
and the size of the gap (hence the spatial size of the 
plateau) increases with magnetic field. One can even 
imagine that for a fixed flux there will appear a sequence 
of FQH states with larger even denominators and thus 
with smaller densities all the way up to the MI-SF phase 
boundary, however their size will be much smaller and 
they may not be discernible at all. Other probes for the 
FQH states may be noise correlations in time-of-flight ex- 
periments, measurement of the Hall conductance for the 
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mass current in a tilted lattice, or Bragg spectroscopy 

§an : 

A major impediment to observing these states is the 
need to reduce the temperature to below the scale of the 
gap, which is a fraction of the hopping matrix element t. 
Such temperatures are currently hard to reach reliably. 

VI. SUMMARY 

In summary, we have predicted that experiments on 
bosons in rotating lattices (or in lattices with an artificial 
gauge field) will see a phase where the excitations on top 
of a Mott insulator form a bosonic fractional quantum 
Hall state. We base our prediction on a set of variational 
calculations, supplemented by exact diagonalization of 
small systems. We find that the MI + v = 1/2 Laughlin 
state has a lower energy than the Gutzwiller mean-field 
vortex lattice when the density of excess particles/holes, 



e = N/L 2 , is chosen appropriately (e = av = a/2), and 
the hopping t is sufficiently small compared to the in- 
teractions U. In this regime we find that the overlap 
between the exact ground state and the Laughlin state 
is as large as 96%, but the overlap with the superfluid 
is smaller than 10%. We produced a phase diagram 
(Fig. Ql, illustrating where this novel phase should be 
found at low temperatures. 
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